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1. 


Introduction 


This  report  represents  the  results  of  the  first  year  of 
a  scouting  program  of  research  in  foundations  of 
mesomechanics . 

Three  main  directions  of  this  program  outlined  in  the 
initial  proposal  were  persued: 

a)  Employment  of  Weyl ' s  geometry  as  a  model  of  material 
space  geometry.  The  potential  of  the  approach  is 
demonstrated  m  a  specific  example:  we  derive  the  classical 
equations  of  thermoelasticity  on  pure  geometrical  ground. 
This  result  is  obtained  by  B.  Kunin  and  A.  Chudnovsky  and  is 
the  content  of  Section  2. 

b)  Development,  of  experimental  techniques  of  studying 
metric  properties  of  the  material  space.  We  report  our 
first  results  in  quantitative  fractography  in  a  search  for 
strength  parameters  revealed  in  fracture  surfaces.  This 
work  consists  of  two  equally  important  parts:  development 
of  experimental  technique  of  observation  and 
characterization  of  fracture  surfaces  for  various  materials 
and  computer  simulation  or  erratic  surface  geometry.  These 
studies  are  being  conducted  by  A.  Kim,  B.  Kunin  and  A. 
Chudnovsky  and  are  reported  in  Section  3. 


c)  Extension  of  the  least  action  principle  from  classical 
mechanics  to  fracture  mechanics.  This  is  the  first  step  of 
application  of  the  Lagrangian  formalism  to  the  continuous 
damage  mechanics.  At  first,  an  introduction  of  a  crack  tip 
Lagrangian  was  thought  of  as  a  training  ground  for 
addressing  the  main  question  of  continuous  damage  evolution. 
However,  as  the  work  advanced,  we  realized  the  importance  of 
the  result  for  understanding  and  description  of  fracture 
phenomena:  a  r.ew  dynamic  crack  propagation  equation 
resulted  from  -his  approach.  In  the  derivation  of  a 
continuous  damage  evolution  equation,  the  dynamic  crack 
propagation  equation  for  an  isolated  crack  will  play  the 
same  role,  as  the  Newton's  equation  plays  in  the  statistical 
mechanics.  The  derivation  of  the  dynamic  crack  growth 
equation  was  done  by  B.  Gommerstadt  (Northeastern 
University,  Boston)  and  A.  Chudnovsky  and  is  presented  in 
Section  4. 

In  Section  5  we  outline  future  work  in  each  of  the 


above  directions. 


2 •  Interpretation  of  Thermoelasticity  as  a  Realization  of 

Weyl 1 s  Geometry. 

The  concept  of  damage  is  central  to  mesomechanics .  In 
general,  damage  may  be  understood  as  a  deviation  of  material 
morphology  from  its  perfect  (ideal)  state.  Thus,  for 
example,  damage  in  a  homogeneous  continuum  can  be  visualized 
as  distributed  heterogeneity  such  as  voids,  microcracks,  new 
phase  inclusions,  etc. 

Usually  a  damage  parameter  is  introduced  as  an  internal 
variable  within  the  framework  of  continuum  mechanics.  This 
approach  is  employed  by  the  majority  of  researchers  in  the 
field  of  damage  mechanics  for  the  last  thirty  years  (see  a 
review  [1]  and  its  update  [2]).  Here  we  explore  an 
alternative  approach  to  modeling  of  damage.  Following 
classical  works  in  dislocation  theory  [3-5],  we  consider  the 
internal  (material)  space  parallel  to  the  external 
(laboratory)  space  and  the  corresponding  geometries.  Metric 
properties  of  the  internal  geometry  are  directly  related  to 
the  material  morphology.  Therefore,  damage,  as  a  change  in 
the  morphology,  results  in  a  change  in  metric  properties  of 
the  material  space.  In  our  exploration  of  geometrical  means 
of  damage  characterization,  we  begin  with  the  Weyl's 
geometry,  since  this  is  a  minimal  deviation  from  the 
Euclidean  geometry  possessing  the  necessary  degrees  of 
freedom.  The  most  lucid  case  of  variation  of  a  metric  in 


the  material  space  is  thermoexpansion.  Thus  we  reconsider 


the  conventional  thermoelasticity  as  a  realization  of  the 
Weyl's  geometry  to  test  the  approach. 

Below  we  define  internal  and  external  metric  fields 
associated  with  a  thermoelastic  solid.  The  internal  metric 
is  intrinsic  to  the  solid  (at  a  given  temperature),  while 
the  external  metric  is  induced  by  an  embedding  of  the  solid 
into  the  laboratory  space.  Strain  as  a  local  measure  of 
deviation  from  a  stress-free  state  is  then  defined  as  a  the 
difference  between  the  two  metrics.  Compatibility  equations 
for  strait,  fellow  then  from  an  apparent  fact  that  the  solid 
with  its  scaied  metric  associated  with  a  non-homogeneous 
thermoexpansion  (a  special  case  of  Weyl’s  Geometry)  is 
situated  in  the  Euclidean  laboratory  space.  The  Hook's  law 
as  a  constitutive  equation  together  with  the  equilibrium 
equations  complete  the  system  of  equations.  This  system  is 
shown  to  be  equivalent  to  the  standard  equations  of 
uncoupled  thermoelastostatics . 

In  what  follows,  p  will  be  material  coordinates 

M 

frozen  within  the  solid,  and  will  be  rectilinear 

Cartesian  coordinates  in  the  laboratory  space.  It  is 

convenient  to  assume  that  5^  were  introduced  as  coinciding 

with  the  laboratory  coordinates  in  the  state  of  a 

homogeneous  temperature  distribution  t  and  the  absence  of 

o 

external  forces  (reference  state).  It  is  assumed  that  a 
thermodynamical  state  of  an  elastic  solid  is  completely 
characterized  by  a  pair  {u,  T),  where  =  x^)  -  £ ±  is 


the  displacement  from  the  reference  state  configuration  and 
T  =  T(|J  is  the  deviation  of  temperature  distribution  from 
the  reference  state  temperature  TQ.  Thus  the  state  {0,0} 
is  the  reference  state. 


To  remain  within  the  scope  of  linear  elasticity,  we 
consider  small  displacement  and  thermoexpansion  only. 


The  external  metric  g  ^  on  the  solid  in  a  state  {u, 
T}  describes  lengths  (and  angles)  of  (and  between) 
infinitesimal  linear  elements  d£_  as  measured  in  the 
laboratory  coorarnate  system  (so  the  external  metric  is 
always  Euclidean!.  Thus  it  is  the  displacement  vector  u 
that  determines  g.  In  fact. 
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neglecting  the  quadratic  in  u  terms.  (Here  a  subscript  i 
after  the  comma  denotes  the  partial  derivative  in  ?i> 
etc.).  By  definition, 


e .  . 

Then 


(2) 


+  2e  .  . 
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(3) 


The  internal  metric  G  .  on  the  solid  in  the  state 

ij 

(u,  T}  describes  lengths  (and  angles)  of  (and  between) 
infinitesimal  linear  elements  d£  as  measured  at  a  point  £ 
at  the  temperature  T0  +  T  (5)  in  the  stress-free  state. 
The  latter,  for  a  general  temperature  distribution,  cannot 
be  achieved  even  if  the  external  forces  are  turned  off. 

The  correct  procedure  is  therefore  to  cut  the  solid  into 
sufficiently  small  pieces,  for  each  of  which  the  conditions 
can  be  met  (up  to  infinitesimals  of  higher  orders).  Since 
in  the  vicinity  of  each  point  the  temperature  distribution 
is  homogeneous  and  homogeneous  increase  of  temperature 
causes  isotropic  expansion  with  the  factor  1+  aT,  then, 
neglecting  the  terms  quadratic  in  aT, 

Gij  *  <1+2“T,4iJ  (4) 


(Kronecker's  delta  represents  there  the  internal  metric  in 
the  reference  state  due  to  the  assumed  choice  of  the 
coordinates  5;  ) .  The  vector  field  V  needed  together  with 


the  metric  G  for  defining  a  Weyl's  geometry  is  equal  to 
the  gradient  of  the  temperature  field:  V  =  T,  . 

The  strain  tensor  in  the  state  {u,  T}  is  defined  as 


e . . 
ij 


e  .  . 
ij 


aT<S  .  . 
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(5) 


From  the  definitions  of  g  and  G..,  it  follows 

ij  iJ 

that  stresses  are  related  to  e  only:  a  =  n(t)  and, 

moreover,  7 (_o )  =  o.  Assuming  The  Hook's  Law  approximation 
as  well,  as  independency  of  the  elastic  constants  C  on 
temperature,  we  can  write 


a  .  . 
V] 


Ekl 


The  equilibrium  equations  are 


a  .  .  .  =  f  . 

1  J  *  J 

where  f  is  body  force. 


(6) 

(7) 


The  missing  compatibility  equations  appear  after  we 
compute  the  curvature  tensor  R  for  G  in  two  ways  (eqn. 
( 11 )  below) . 


On  one  hand,  we  substitute  G  by  e  -  2z  (from 

ij  ij  ij 

(5))  and  neglecting  the  terms  of  higher  than  the  first  order 
in  e  (since  e  is  small),  we  arrive  at 


Jir 


£  i  1  >  kr  £ik,jr 
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jr,ik 


kr ,  ij 


(8) 


f that  g  disappeared  is  related  to  the  fact 
a  Euclidean  metric). 


g 

ij 


is 


On  the  other  hand,  we  substitute  G  .  by  (1  +  2aT)6i, 

ij  1J 

(from  (4))  and  neglecting  the  higher  order  terms  in  aT >  we 


R  =  -  x  (  5  .  .  T  -  5 ..  T  .  -  6  .  T  ..  +  «.  T,  ,)  .  (9) 

kj  ir  ij  ik  >jr  jr  ,  lk  kr  ij 

Comparing  (9)  to  (8)  we  get  the  compatibility  equations  for 
strain,  though  written  in  an  inconvenient  form: 


ij  ,kr  ~  ’  Lk, j  r  " jr, ik  kr,ij 


-  5  T  -  6  T.  +  6,  T  .  . )  . 
' ’ ' i i ‘ , kr  ik  ,jr  jr  ik  kr  , ij 


The  expressions  above  are  antisymmetric  within  the 

pairs  of  indices  kj  and  ir,  therefore  contracting  with 

e  e  results  in  an  equivalent  equation  with  only  two 

lkj  rair 

free  indices 


4Rotln  -  =  4a(T,kk1 


-  T,  ) 
m  lm 


or,  since  we  assume  the  temperature  distribution  T 
stationary  and  thus  T,kk  =  0, 


Rot£ui  ~  ~  lT7m  * 


Here  the  operator  Rot  is  defined  by 


KUtj  p  t.  .  O  ..  .  UJI 

cm  ~  tkj  mir  j  r , ki 

*  is  the  completely  antisymmetric  tensor  with  the 

components  in  Cartesian  coordinates  equal  to  +1  and  -1 
for  even  and  odd  permutations  of  1,2,3  respectively  and  zero 
for  ail  other  combinations  of  indices. 


for  any  tensor  p 

ii 

Finally,  the  complete  system  of  equations  is  ((7),  (6), 

( 11)  ) 


together  with  coundary  conditions,  which  are  assumed  here  to 
be  in  terms  of  forces.  If  interested  in  the  displacement, 
one  comhines  (2)  and  (5) 


and  solves  (17)  for  u  (upon  finding  £  from  (14-16). 

For  an  isotropic  material,  C  .  =  2y6.,5.„  + 

ijkc  ik  jZ 

^  „  and  it  is  easy  to  see  that 

ij  kX 

equations  (14)  -  (17)  are  equivalent  to  the  conventional 

system  of  equations  of  the  uncoupled  quasistatic 
thermoelasticity.  Indeed,  the  conventional  system  is  (Eqs. 
(2.6.2)  -  (2.6.4)  in  [6]) 
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(18) 
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f . 
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cr..  =  A6..c,,  +  _^e..  -  13A  +  2jj  )  5..aT  (19) 
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e  .  . 
ij 


2  (ui,i 


(20) 


Here  the  displacement  field  u  is  considered  a  primary 
variable  (  y  is  a  derived  variable).  Then  e  identically 
satisfies  (16)  (with  (2),  (5)  taken  into  account ) .  In  the 
preceeding  exposition,  the  strain  £  is  the  primary 
variable  and  - r.e  displacement  can  be  reconstructed  using 
(  17)  . 


The  geometry  utilized  above  is  a  very  special  case  of 
Weyl's  geometry.  General  Weyl's  geometry  is  more  rich  and 
appears  adequate  for  modelling  distributed  damage. 
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Fracture  Surface  Characterization 


-5  . 

It  is  generally  understood  that  material  toughness  and 
specific  fracture  energy  are  higher  for  more  erratic 
fracture  surfaces.  Fractal  geometry  suggests  parameters  to 
quantify  erraticism  of  crack  surfaces.  One  of  them  is  the 
crack  diffusion  coefficient  used  in  the  'crack  diffusion 
model'  of  brittle  fracture  [1],  Other  characteristics  are 
supplied  by  fractal  studies:  fractal  dimension  D  of  a 

fracture  surface  'profile)  together  with  basic  patterns 
l "generators" )  on  various  scales. 

To  study  the  correlation  between  D  and  material 
toughness,  one  needs 

a)  to  determine  which  of  the  existing  techniques  of 

experimental  evaluation  of  D  is  most  suitable  for  a 

given 

material , 

b)  to  analyze  the  consistency  of  an  experimental  procedure 
under  consideration . 


1  • 1  Some  Methods  of  Finding  Fractal  Dimension 

There  are  two  approaches  to  finding  the  fractal 
dimension  ( D )  of  a  fracture  surface:  Profile  Analysis  and 
Surface  Analysis.  Figures  1  and  2  represent  the  two 
approaches  schematically. 


i 


3 • 2  Outline  of  Performed  Studies 


Profile  Analysis 


Richardson's  method  uses  a  sequence  of  measurements  of  the  "length”  of 
a  profile  utilising  a  compass  of  decreasing  opening,  i.e.  a  sequence 
of  lengths  of  increasingly  fine  polygonal  approximations  of  the  profile 
Also  known  as  the  "measuring  stick  technique”. 

Minkowski's  method  of  determining  D  of  a  planar  profile  utilizes 
a  sequence  of  areas  of  thinner  and  thinner  strips,  each  obtained  as 
a  result  of  uniform  thickening  of  the  profile.  The  ratios  of  the  areas 
to  the  strip  widths  are  analogs  of  polygons  in  the  Richardson's  method. 
The  Fourier  Transform  method  utilizes  the  sequence  of  Fourier  coef¬ 
ficients  of  the  profile,  specifically,  the  rate  of  their  convergence 
to  zero.  The  lengths  of  the  graphs  of  harmonics  in  the  profile's 
Fourier  series  are  analogs  of  the  lengths  of  polygons  in  the  Richard¬ 
son's  method. 


Figure  1 


Surface  Analysis 


"Slit  island  analysis"  of  Mandelbrot  and  Passoia  (thick  arrows  above) 
utilises  a  D-dependent  relation  between  areas  and  perimeters  of  slit 
islands,  as  the  latter  appear  and  grow  in  the  plane  of  polish. 

Minkowski's  method  of  determining  D  of  a  surface  in  3-space  uses  a  sequence 
of  thinner  and  thinner  layers,  each  obtained  as  a  result  of  a  uniform 
thickening  of  the  surface.  (The  ratio  of  the  volume  of  the  layer  to 
its  thickness  is  an  approximation  to  the  surface's  "area".) 


Figure  2 


The  experimental  program  conducted  in  the  first  year 
consisted  of  a  few  steps. 

a)  Analysis  of  which  of  the  existing  experimental 
techniques  is  most  suitable  for  a  particular  material  (see 
Figs.  1,  2).  Our  model  materials  were:  graphite  (which 
exhibits  typical  features  of  fracture  surfaces  associated 
with  brittle  failure  of  polycrystalline  materials)  and 
polyethylene  (with  highly  fibrillated  fracture  surfaces 
typical  for  many  engineering  plastics). 

b)  Application  of  particular  techniques  to  the  model 
materials.  The  "slit  island"  method  was  applied  to  graphite 
and  the  Richardson's  plot  method  to  both  graphite  and 
polyethylene . 

c)  Computer  simulation  of  fractal  trajectories  by  means  of 
Fourier  synthesis.  Comparison  of  the  simulated  trajectories 
with  the  observed  profiles. 

In  Figs.  1,2  we  summarized  most  popular  techniques  of 
fractal  dimension  evaluation.  There  are  two  techniques 
applied  directly  to  surfaces.  The  absorption  method  leaves, 
in  our  opinion,  an  uncertainty  due  to  the  limitations  in 
controlling  the  uniformity  of  coating  thickness.  The  "slit 
island"  method  appears  straightforward  and  has  been 
employed.  It  remains  to  be  seen  how  sensitive  it  is  to  the 
orientation  of  the  plane  of  polish. 


Out 


"tree  methods  applicable  to  fracture 
profiles,  the  Furrier  -Technique  seems  more  promising,  hut  it 
requires  sorr.vcar  ~r.at  v;e  do  not  have  at  present.  The  ether 
two  methods  -.re  'imiiar  in  their  limitations  and  advantages; 
we  employed  -.to  *  er.ardson '  s  plot  method. 

Below  we  report  our  results  in  application  of  the  "slit 
island"  ano  ~  -  to amen1  s  plot  techniques. 


ncrue 


-he  mur;:;  rurrace  was  obtained  from  a  compact 
tension  rrapr.ite  -oecimen.  Figure  3  shows  a  general  view  of 


Fig .  3 

tne  fracture  surface.  To  proceed  with  comparative  studies 
of  different  tecnmques,  multiple  wax  replicas  of  the 
fracture  surface  were  made  and  "shaving"  was  performed  on 


Figure  6  shcwo  -  r.o  Image  Analyzer  screen  with  an  island. 


The  image  c 
( required  by 


digitized  and  its  area  and  perimeter 
ccedure)  are  evaluated  automatically. 


Fis? .  6 

Figure  7  shews  tr.e  observed  area  vs.  perimeter  relationship 


[ dots )  and  ;:s 


ximation  by  a  straight  line. 


Fractal  Dimension 
D-1.2? 


log  P 


The  line's  slope  =  1.57  yields  the  fractal  dimension  of 
the  fracture  surrace:  D  =  1  +  2/k  =  2.27. 

-  ■  Richarcscr. ' s  Plot  Technique 

For  grapnitc ,  the  fracture  profiles  were  obtained  from 
tne  specimen  sr.cwr.  .. r.  Fig.  3  above.  Wax  replicas  were  used 
again  to  prcsorv'  -.no  specimen.  A  typical  profile  is  shown 
in  Fig .  .3  . 


Fig.  8 

Figure  9  shows  the  observed  profile's  "length"  L  as  a 
function  of  the  compass  opening  (dots)  together  with  the 


straight  line  approximation  of  the  relation. 


a  2  3  4  s 

log  E 


F  i  s> .  9 

The  j.cpe  •:  =  -0.23  of  the  straight  line  yields  the 
fractal  dimers  ter.  if  the  profile:  D  =  1  -  k  =  1.23,  i.e.  the 
dimension  l-d  =  2.12  of  the  fracture  surface. 


M 


For  t  ratioue  fracture  surface  in  polyethylene  (Fig. 
10),  fiaures  ll, 12  show  a  typical  fracture  profile  (parallel 
to  the  crac.-:  front)  ana  the  prof  iie-length-vs .  -compass- 
opening  relatrcnsmp . 


The  slope  k  =  -0.08  of  the  straight  line  in  Fig.  12  yields 
D  =  1.08  for  the  profile. 

3.5  Profile  Simulation 

To  simulate  the  profile  from  Fig.  8  the  following 
relation  between  the  fractal  dimension  of  a  curve  and  the 
rate  of  decrease  of  the  curve's  Fourier  coefficients  [5]. 
If  a  funct  icr.  f 1  x )  is  expanded  in  a  Fourier  series 

~r 

f  ( x )  =  a  *  )  a  sm(kx  +  0) 

0  k  k 

and  a  ~  k"1  as  k^  «  ,  then  the  graph  of  the  function 

k 

f(x)  has  ( generically )  the  fractal  dimension  D  =  2-a. 

We  generated  sequences  of  coefficients  a^  by  calling 

a  random  number  generator  (for  each  k),  which  supplied  us 

with  a  random  number  r^  equidistributed  between  0  and 

1,  and  then  put  a  =  ( 2r  -  l)k_ct  ,  so  that  a  is 

k  k  K 

“CX 

equidistributed  between  ( -k  )  and  k 

Figure  13  shows  a  curve  obtained  for  D  =  1.23  (a=  2-D  = 

.77).  The  Fourier  series  was  truncated  to  500  terms;  a 
v  ° 

and  all  0  were  put  equal  to  zero  for  simplicity, 
k 


i  C  .  13 

*  .  6  T cr.c*  !  r  ~ 

Ext:r.c:vo  :<per  imentation  leads  us  to  the  following 
conclusions . 

We  possess  -.uocuate  techniques  of  evaluation  of  fractal 
dimension . 

2.  Fractal  dimension  alone  is  insufficient  for  adequate 
simulation  ef  fracture  profiles  (Fig.  13  clearly  does  not 
mimic  the  anginal  profile).  One  needs  additional 

characteristics:  casic  patterns  ("generators"),  scale 
dependency,  etc. 
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4.  Dynamic  Crack  Growth  Equation 
A.l  Introduction 

The  kinetics  of  fast  crack  propagation  has  been  discussed  since  the  pio~ 

? 

neering  work  by  Mott*  for  the  last  40  years. 

A  choice  of  a  fundamental  assumption  which  leads  to  dynamic  crack  propa- 

2,3 

gation  equation  is  the  Key  point  of  the  discussions.  Proposed  by  Eshelby 
and  presently  widely  accepted  way  to  derive  dynamic  crack  growth  equation  "is 
to  suppose,  along  with  A.  A.  Griffith,  that  the  energy  released  at  the  top  is 
all  used  to  provide  'he  surface  energy  (y)  of  the  faces  of  the  crack.  This 
gives  the  equation  ut  motion"  (quoted  from  ref.  2): 

G?(£,c*'  -  2-,  .  (1) 

d 

where  G^  is  the  dynamic  energy  release  rate  with  respect  to  the  crack  tip 

advance  and  £  and  £  are  the  crack  length  and  velocity  respectively.  (2~D 

problem  is  considered).  Achievements,  limitations  and  difficulties  of  dynamic 

crack  growth  theories  as  well  as  the  comparison  with  experimental  data  are 

4 

quite  adequately  presented  in  several  reviews,  for  instance  Erdogan  and 

n  5 
Rose  . 

Here  we  emphasize  one  problem  which  is  associated  with  the  nature  of 
dynamic  crack  extension.  It  follows  from  elastodynamic  solution,  that  the 

dynamic  energy  release  rate  does  not  depend  on  crack  tip  accelerat i on.  Thus, 
if  the  crack  tip  is  regarded  as  a  particle,  then  it  is  the  one  which  exhibits 
no  inertia.  It  is  a  puzzling  result.  Indeed,  a  crack  can  be  formally  repre- 
sented  as  a  continuous  array  of  dislocations  moving,  when  the  crack  extends. 
It  is  also  well  known  that  every  individual  dislocation  possess  an  "inertia". 


However  a  pile  of  them  loses  this  feature.  Esheiby  suggested  a  possible 


explanation:  "In  a  limit  of  a  large  number  of  dislocations,  their  mutual 

interactions  entirely  dominate  their  individual  properties".  This  guess  still 

remains  unconfirmed.  Beside  this  theoretical  difficulties,  there  are  experi- 

7  8 

mental  evidences,  that  a  dynamic  crack  exhibit  inertial  behavior  .  One  of 
such  observations  stimulated  present  study.  We  consider  an  alternative 
approach  to  the  derivation  of  dynamic  crack  growth  equation:  we  define  a  crack 
tip  Lagrangian  employing  e 1  a s t odynam i c  solution  and  then  apply  the  Least 
Action  principle  to  select  the  actual  crack  tip  motion  from  an  ensemb 1 e  of  all 
pos  s i b  1  e  ones  . 


4.2  Kinetics  of  Crack  Propagation 

We  limit  our  consideration  to  a  rectilinear  crack  propagating  through  a 
homogeneous,  isotropic,  i  inear  elastic  solid. 

Let  us  consider,  the  conventional  Lagrangian  density  of  el  as t odynami cs : 


L 


f 


(1) 


where  k  =  y  p'x  lu.u.  -  stands  for  kinetic  energy  density  f  ^  Ui  j  ^ijk£ 


u  ,  is  the  strain  energy  density. 


x  ,  u.,  u.  .  and  u.  are  the 
o’  1  i.j  i 


P  and  C  are  the  density  and  the  elastic 
moduli  tensor,  x.  is  Cartesian  coordinate, 

i 

time,  displacement  vector,  displacement  gradient,  and  the  velocity  components 
respectively.  Latin  and  Greek  subscripts  take  values  i,j,k,  =  1,2,3  and  or, 3 


=  0, 1.2,3. 

The  global  Lagrangian  of  an  elastic  solid  V  can  be  obtained  by  inte¬ 
grating  over  the  volume  V; 


Then  conventional  action  results  from  the  integration  of  L  with  respect  to 


time. 

Variation  of  the  action  with  respect  to  displacement  leads  to  conven- 

The  following  balance 

J  \dUi,i/ 

equations  also  hold: 

bt\  AP.wJ  (3) 


tional  equations  of  e  1  as  t  odynami  cs  pu.  =  i)  .  (  -  — _ 

1  J  3u.  . 

\  i.J 


ftx  j  exp 


ax 


Here  (Vbx) 


^./exp  means  the  derivation  with  respect  to  x  when  the  rest  vari“ 

or 

ables  are  held  fixed  and  P  is  the  4X4  e 1  as t odynami c  energy  momentum  ten¬ 


sor,  introduced  by  h shelby 


1.6 


crll 


■>u 


;< ,  d 


u  < > ,  u 


The  components  of  energy  momentum  tensor  can  be  expressed  in  terms  of  kinetic 
and  strain  energy  densities,  stress,  strain  and  displacement: 


p.  - 
1 J 

(  f  -  k )  i'i 

p, 

Li  -  IJ  . 

k  o 

i  i  ,k 

P  .  = 

-  a  .  .  u  . 

°J 

i  J  i 

P 

00 

f  +  k 

i  j  k  ,  j  k  ,  i 


(4) 


We  consider  a  crack  in  a  plane  x^.x^  w i t h  a  straight  front  along  parallel  to 
*3  moving  in  x^  direction.  Then  the  crack  front  location  is  determined  by 
x j-coordi nat e  x^=  £(x^),  where  t  is  the  crack  length,  xq  stands  for  time. 
We  describe  the  crack  as  a  singularity  (zero-values  in  density  and  elastic 
moduli  along  the  crack)  Thus, 


(5) 


p  (x  )  =  p (x . -£) 
a  1 

C (x  )  =  C  (  x  .  ~£) 
-  u  —  L 


and  consequently 


From  eqs 


dl 

dx 

H. 

1 

exp 

(2)  . 

(3) 

and  (0)  we  arrive  at 

r 

b  L 

Id 

= 

dv 

->x 

J 

V 


Applying  it  we 


introduce  a  crack  Lagrangian 


LC  (1,1): 


(6) 


(7) 


l 


r  'ic  c  . 

j  ^  <*  '  L  <*> 


(8) 


where 


o£c 

-H 


i  m 

v-<0 


;,rfpiadv  and 


(9) 


£C  (1)  can  be  evaluated  as  (Z,Z)  at  Z-* 0 

The  definition  (9)  is  based  on  the  assumption  that  the  elastic  field 
within  an  i nf i n i t es i na l  vicinity  of  the  crack  tip  only  governs  crack  propaga¬ 
tion.  Apparently,  L°  is  only  a  part  of  the  total  Lagrangian  of  a  solid  with 

the  crack.  Now  we  apply  the  Lagrangian  technique  to  formulate  the  dynamic 

c 

crack  motion  equation  as  an  Eul er_Lagrange  eq.  associated  with  L  : 

d  ;>Lc  ->Lc 


dX 


■)Z 


=  ~2y. 


(10) 


0  bZ 

where  the  specific  energy  of  new  surface  formation  2y  appears  as  the  energy 


sink. 


To  specify  eq.  (10)  we  need  to  account  for  the  elastic  field  near  the 


crack  front.  !n  this  region  the  asymptotic  displacement  field  can  be 
expressed  as  u(X^.  X0,  X  )  =  u(X^  -  x^,  X^^and  consequently 


It  can  be  shown  (see  for  instance  [5])  that  the  integrand  of  eq.  (9)  can  be 
presented  in  the  divergence  form: 


'P 


■)P 
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W, 


IX 


*=  ^  H 
J  1 J 


(11) 


where  H.  .  -  i[*k)  u  .u,  ..  Then  following  the  conventional  notations 

i  j  i  j  k  j  k  ,  i  6 

we  introduce  a  dynamic  energy  release  rate^. 


G1  <£,£>  ~  ll(n 


v  0  ' 


lHi j  "j 

Ov 


(12) 


where  is  the  unit  outer  normal  to  the  boundary  SV 

3 

As  has  been  shown  can  be  factorized: 


Gj  %V  =  HV  gl<£) 


(13) 


where  G(£)  is  the  conventional  static  energy  release  rate  per  unit  crack 
advance,  g(£l  is  the  velocity  factor. 

Upon  substitution  (13)  into  (8)  we  obtain: 

L  C  (£,  £  )  -  g  (£)  n  (£)  .  (1A) 

where  (£ )  is  the  potential  energy  of  the  solid  associated  with  the  crack. 

c 

The  term  l  on  (8)  vanishes  since  l°^£,£^  tends  to  zero  with  £-»0  •  The 

latter  follows  from  Kostrov's  solution  for  dynamic  stress  intensity  factor  for 

9 

a  crack  suddenly  appearing  in  an  initially  continuous  solid. 


Making  use  of  CIA)  and  (10)  we  arrive  at  the  following  crack  propagation 
equat i on 

m  (l, hi  =  F  (1,1)  ,  05) 

here  we  designate  'm'  as  a  crack  inertia  coefficient  (CIC) 

m (lfb  -  < t>(l)  n  (t)  ,  06) 

and  the  crack  driving  force 

F  (lfb  -  <t>lb  C  l£)  -  2y(h  07) 

.  1  d2e(M) 

$>(£)  =  - — 

c  dM~ 

dg  (M) 

=  s(M)  ‘  M  “dT 

M  =  i/c 

4.3  Discussion  and  Conclusions 

•  • 

Since,  m  and  F  do  not  depend  on  Z9  then,  according  to  eq.  (15)  a 
crack  propagation  appears  as  a  pseudoparticle  motion  with  a  "mass"  being  pro¬ 
portional  to  the  potential  energy  associated  with  a  crack.  The  crack  driving 
force  represents  the  excess  of  energy  release  rate  above  the  required  fracture 
energy  2y . 

The  expression  C^(£)  for  various  loading  conditions  and  specimen  config- 

« 

urations  is  readily  available  [10].  The  velocity  factors  g(£  )  and  conse- 
quently  ^C£)  and  <t>(£)  require  a  solution  of  elas todynamic  problem  and 
have  an  exact  solution  for  only  a  few  specific  loading  conditions. 


Ail  illustrative  example  of  the  dynamic  anti-plane  crack  growth  equation 
is  presently  in  preparation. 

We  have  two  concluding  remarks  to  make. 

1.  The  crack  propagation  equation  needs  an  experimental  examination. 

2.  It  is  well  known  that  fracture  surfaces  created  at  different  crack  speeds 
exhibit  different  roughness.  In  a  continuum  based  model  outlined  in  this 
section,  the  roughness  of  a  fracture  surface  can  be  accounted  for  by  a 
dependency  of  the  specific  fracture  energy  on  the  crack  speed:  y  =  7(£)  • 
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5 .  Future  Work 


V 


We  consider  the  following  tasks  as  a  logical 
continuation  of  the  research  outlined  above. 

a)  To  formulate  the  kinematics  of  a  continuum  theory  of 
defects  (state  variables,  compatibility)  on  the  basis  of  the 
Wey 1 ' s  geometry . 

b)  To  develop  a  methodology  of  experimental  evaluation  of 
the  parameters  of  the  Weyl's  geometry  (fractal  type 
characteristics  at  various  scales) . 

c)  To  expand  the  Lagrangian  formalism  to  media  with 
evolving  damage;  to  combine  the  resulting  constitutive  laws 
with  the  obtained  compatibility  equations  and  the 
equilibrium  equations  and  examine  the  entire  model 
experimentally . 

d)  To  develop  and  examine  experimentally  a  model  for  defect 
nucleation  to  account  for  processes  occurring  in  the 
material  before  the  Lagrangian  formalism  takes  over. 


